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We present a comparative analysis of a set of GaAs-based THz quantum cascade lasers, based on longitudinal- 
optical phonon scattering depopulation, by using an ensemble Monte Carlo simulation, including both carrier- 
carrier and carrier-phonon scattering. The simulation shows that the parasitic injection into the states below 
the upper laser level limits the injection efficiency and thus the device performance at the lasing threshold. 
Additional detrimental effects playing an important role are identified. The simulation results are in reasonable 
agreement with the experimental findings. 
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Quantum cascade lasers (QCLs) have an enormous po- 
tential as compact and efficient THz sources, promising 
mW-level output power. Recently a new scheme based 
on longitudinal-optical (LO) phonon scattering depop- 
ulation has been successfully demonstrated^ - — designed 
for high operating temperatures and reaching up to 164 K 
in pulsed and 117 K in continuous wave (cw) operation. 4 
In these structures, the collector state is separated from 
the next lower level by at least the LO phonon energy, 
enabling efficient depletion of the lower radiative state by 
LO phonon scattering and reducing thermal backfilling. 
This overcomes the limitations associated with the con- 
ventional bound-to-continuum and chirped superlattice 
designs, where thermal backfilling of the lower radiative 
state prevents high temperature performance. 

In this letter, we theoretically investigate three 
GaAs/Al .i5Gao.85As THz QCL structures based on the 
above scheme and fabricated with high-confinement low- 
loss double metal waveguides. An experimental compar- 
ison between these structures was already carried out^ 
and motivates our theoretical comparative investigation. 
For the analysis, an ensemble Monte Carlo (MC) method 
is used,— ~— which has successfully been applied to the the- 
oretical investigation of such structures^— The subband 
energies and wavefunctions of the structure are calcu- 
lated with a Schrodinger-Poisson solver, coupled to the 
MC simulation. The essential scattering mechanisms due 
to electron-electron, electron-LO phonon, and electron- 
acoustic phonon interactions are considered as described 
in Ref. 6. Intercarrier scattering is implemented based 
on the Born approximation; 11 ' 12 also taking into account 
screening effects^ In the simulated structures, intercar- 
rier scattering plays a major role for the intrasubband dy- 
namics, while carrier transport is governed by LO phonon 
scattering for the temperature and doping level consid- 
ered. Non-equilibrium phonon distributions are explicitly 
taken into account, and the Pauli's exclusion principle is 
also included. Due to the relatively low electron densities, 
the influence of such mechanisms is negligible. Periodic 
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FIG. 1. (Color online) (a) Conduction band profile and 
probability densities of a 2.8 THz QCL based on a phonon 
depopulation scheme at 11.3kV/cm. The layer thick- 
nesses of one period, starting from the injection barrier, 
are 56/93/36/160/39/67/25/81 [A]. The underlined layer is 
doped at n = 1.9 x 10 16 ci 
period is n,2D = 3x 10 10 cm 
of the applied field for two periods of the 2.8 THz structure 
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boundary conditions are used, i.e., electrons coming out 
from one side of the device are automatically injected 
into the equivalent level on the opposite sideJ Fig. 1(a) 
shows one period of the conduction band structure of 
the 2.8 THz QCL£ calculated at 11.3kV/cm; the band 
profiles of the other structures discussed here are simi- 
lar. The upper and lower levels of the lasing transition 
around 2.8 THz are labeled 5 and 4, respectively. Effi- 
cient depopulation of level 4 is achieved via the collec- 
tor state 3, which is in turn separated from the levels 2 
and 1 by roughly the optical phonon energy (36meV), 
thus enabling efficient depletion to these levels by op- 
tical phonon scattering. Parasitic channels counteract- 
ing the onset of population inversion between levels 5 
and 4 are provided by the injection from levels 1' and 
2' (i.e., the replica of levels 1 and 2 situated in the ad- 
jacent period on the right side of Fig. 1) into levels 1 
- 4, and by scattering processes from the upper laser 
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FIG. 2. (Color online) Simulation results for the 2.8 THz 
structure as a function of the applied field (E). The crosses 
mark the sampling points; the lines are guide to the eye. 
(a) Relative population difference (solid line) and inversion 
(dashed line) between the upper and lower laser level. The 
circles mark the positions of the anticrossings 5-2' and 5 - 
1'. (b) Current density, plotted as a function of E. 

level to the states below. In Fig. 1(b), the level ener- 
gies of the 2.8 THz design are shown as a function of the 
applied electric field. The strong anticrossings between 
level 5 and the injection levels 2' and 1' around 9.9 and 
llkV/cm facilitate a selective injection into the upper 
laser level, leading to a reduction of the parasitic chan- 
nels. Around 7.6kV/cm, there is an additional anticross- 
ing between levels 1' and 4, causing a strong parasitic in- 
jection into the lower laser levels The energy difference 
between the collector state 3 and the injector state 2(1) 
increases from 28meV (36meV) for 7kV/cm to 40meV 
(44meV) for 13kV/cm, enabling an efficient LO-phonon- 
induced depletion of the lower laser level especially for 
biases > lOkV/cm. Figure 2 shows MC simulation re- 
sults for the 2.8 THz structure at a lattice temperature 
of 100 K. The sampling electric fields were chosen as to 
avoid the excessive current spikes emerging near narrow 
anticrossings in MC simulations^ It should be pointed 
out that sharp anticrossing features are an artifact of the 
Schrodinger solution in the absence of broadening mech- 
anisms (such as those provided by surface roughness) . If 
a more precise quantum mechanical approach were used, 
based for example on density matrices^— or nonequilib- 
rium Green's funcionsii, much smoother features would 
be observed, that would nevertheless not drastically mod- 
ify our analysis. In Fig. 2(a), the difference and ratio of 
the relative population in the upper (p§) and lower (^4) 
laser levels is displayed, with A7554 = p$ — p±. The cur- 
rent density in Fig. 2(b) exhibits a parasitic peak (not 
shown) around the 1' - 4 anticrossing and a local mini- 
mum at 8.25kV/cm. It reaches its maximum around the 
two anticrossings of level 5 with the injector states V and 
2', where the injection into subband 5 is most efficient. 
In addition, the energy difference between the collector 
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FIG. 3. (Color online) Population difference between upper 
and lower laser level as a function of the applied field. The 
circles mark the positions of the anticrossings 5-2' and 5 - 
1'. The inset shows Aps4 as a function of the current density 
close to the onset of inversion. The crosses mark the sampling 
points; the lines are guide to the eye. 

state 3 and levels 2,1 is designed to reach the LO phonon 
energy in this bias range and thus become favorable for 
LO phonon-induced depletion of level 4. These two ef- 
fects combine to yield a large population difference Aj>54 
and thus a high material gain between 10 and 11 kV/cm. 

The calculated inversion peak p^/pi = 5.1 lies above 
the experimentally measured maximum p^/p^ — 2.58 ± 
0.52f2> and also the experimentally measured current den- 
sity peak^ is overestimated by about 40%. This is a 
well known phenomenon for MC simulations of structures 
dominated by LO phonon-induced carrier transport 
and may partly be due to the absence of broadening 
mechanisms in the algorithm, as discussed above. To 
some extent, the deviation can also be ascribed to a 
lower than specified free-carrier density in the sample. 
In addition, the experimental effective current density, 
which has been calculated based on an effective area, 
is smaller than the on-axis value in the gain medium. 
The simulated inversion peak is shifted by — 2.4kV/cm 
relative to the experimental maximum, which occurs at 
11.9kV/cnifi2 and the current density peak is shifted by 
a similar amount. Possible explanations are an additional 
parasitic resistance in the experimental structure or de- 
viations in the growth process. In Fig. 3, the population 
difference AP54 is plotted as a function of the electrical 
field for three THz QCLs operating at 2.8 THz 3 , 3.2 THz 2 
and 3.8THai. The maximum of AP54 in Fig. 3, which 
occurs at an especially high bias for the 3.8 THz structure 
in agreement with experiment ji is for all designs located 
around the two anticrossings of level 5 with the injec- 
tor states. The inset shows Ap^ close to the onset of 
inversion as a function of the current density, for values 
ranging between the minimum value after the parasitic 1' 
- 4 anticrossing peak and 0.4kA/cm 2 . This corresponds 
to an electrical field range of 8.25 — 9.4, 9.15 — 10.4 and 
9.85 - 10.5kV/cm for the 2.8, 3.2 and 3.8 THz design, 
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respectively. In the 2.8 and 3.2 THz structures, the on- 
set of inversion occurs at much smaller current densities 
than for the 3.8 THz one, which is an indication of the 
weaker parasitic channel 1' - 4. These findings are in 
agreement with the experimental results, which yield a 
significantly reduced threshold current for the 2.8 and 
3.2 THz structures as compared to the 3.8 THz design^ 
To further investigate the structures, we introduce the 



TABLE I. Overview over simulation results for the 
2.8/3.2/3.8 THz structure close to the onset of inversion. 



Quantity 

P5 (%) 
P4 (%) 
Vinj,5 (%) 
r?inj,4 (%) 
T]inj,3 (%) 
V^nj,2 (%) 
(%) 

r 5 (ps) 
r 5 ^ 4 (ps) 

T4 (ps) 



J = 0.3kA/cm 2 
7.5/8.2/3.5 
1.6/1.5/1.2 
43.3/42.8/26.1 
14.6/15.2/22.6 
20.9/23.6/17.8 
3.5/3.2/12.2 
14.2/10.7/19.2 
2.75/3.05/2.02 
11.1/11.7/13.7 
0.97/0.88/0.68 
0.65/0.53/0.78 



J = 0.4kA/cm 2 
10.6/13.8/7.9 
2.1/2.5/1.3 
42.9/48.4/40.1 
16.3/18.9/13.6 
23.0/21.9/26.1 
3.4/2.2/9.3 
13.0/7.7/9.9 
2.96/3.43/2.18 
10.1/9.7/11.9 
0.88/0.83/0.71 
0.52/0.49/0.72 



injection efficiency Tfi n j i from a given period into a level 
i of the next-lower period, defined as the ratio between 
the carriers injected into level i and the total amount of 
carriers injected into that period. An overview of the 
simulation results is given in Table I. The low inversion 
in the 3.8 THz structure at J = 0.3 kA/cm 2 can be partly 
attributed to the poor injection efficiency into the upper 
laser level Jftry.S) which is only 26% for the 3.8 THz de- 
sign, as compared to over 40% for the other structures. 
This is due to the strong parasitic injection into the states 
below the upper laser level, as indicated by the increased 
Vinj,i.A- From Table I we can see that also for the 2.8 
and 3.2 THz designs, the parasitic injection channels to 
levels 1-4 play a significant role in the threshold re- 
gion. To obtain a high inversion, long lifetimes T5 and 
T5^4 are desirable^ as well as a small T4 ensuring a 
fast depopulation of the lower laser level. Furthermore, 
the gain is directly proportional to the reduced oscilla- 
tor strength f 5 4 4 T5 in the 3.8 THz structure is reduced 
mainly due to an increased leakage to levels 3 and 1, 
while the values of T4 and /54 are slightly superior for 
this design. At low current densities, the 2.8 THz struc- 
ture shows a higher value of ^54 than the 3.2 THz laser, 
outweighing the somewhat lower inversion and resulting 
in a slightly better threshold performance, as experimen- 
tally observed 4 

We have also calculated the electron-lattice energy re- 
laxation rate 1 ^ Such quantity reaches its peak in the 
lasing region, where the phonon-assisted injection into 
the upper laser level and depletion of the lower laser 
level are the strongest. For the best of the structures 
(2.8 THz), a maximum 1 of approximately 3.8 ps -1 was 



calculated, in reasonable agreement with experiment, 3 
The other two structures show a lower maximum value 
of r^ 1 (3.3 ps -1 ), also in qualitative agreement with 
the experiment. For some of the subbands, the simu- 
lation yields highly nonequilibrium carrier distributions, 
complicating comparison with experimentally extracted 
temperatures? 3 - However, an investigation of this issue 
goes beyond the framework of this paper and will be ad- 
dressed in a different article. 

In conclusion, we have shown that the MC simulation 
allows us to identify the parasitic processes affecting the 
operation of THz lasers, thus providing a meaningful ex- 
planation of the experimental findings. 
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